A FREE-BOUNDARY PROBLEM FOR THE EVOLUTION 
p-LAPLACIAN EQUATION WITH A COMBUSTION BOUNDARY 

CONDITION 



TUNG TO 



Abstract. We study the existence, uniqueness and regularity of solutions of the equation 
ft ~ Apf — div {\Df\''~^ Df) under over-determined boundary conditions / = and 
\Df\ = 1. We show that if the initial data is concave and Lipschitz with a bounded and 
convex support, then the problem admits a unique solution which exists until it vanishes 
identically. Furthermore, the free-boundary of the support of / is smooth for all positive 
time. 



1. Introduction 



Fix a number p > 2. Given a non-negative function /o on M" with positive set Qq, we 
want to find a non-negative function /(x, t) on M" x (0, T) with positive set fl which solves 
the following problem: 



(P) 



ft = \f mn = {f>0} 

/ = and \Df\ = 1 ondnn{0<t<T} 

limt^o/(x,t) = /o(x) VxGM". 



The operator 

Apf = div{\Df\P-^Df) 
is known as the p-Laplacian. In non-divergent form, it can be written as 



(1.1) 



Ap = iDfr^Af +{p- 2)\Dfr^fi,pjj 



Note that the Einstein summation notation was used in the last term. It can also be written 
as 



(1.2) 



Ap=\Dfr\Af + {p-2)f, 



where f^i, denotes the second derivative of / in the direction of = Df /\Df\. 

In the case p > 2, this operator is nonlinear and degenerate at vanishing points of Df. 
When p = 2, it is just the regular Laplacian. 



2000 Mathematics Subject Classification. Primary: 35R35; Secondary: 35K55, 35K65. 

Key words and phrases. p-Laplacian, free-boundary problem, degenerate equation, combustion, regulariza- 
tion, convex domain. 

1 



2 



TUNG TO 



Due to the over-determined boundary conditions / = and \Df\ = 1, the time-section 

= {x G M" I /(x,t) > 0} 

will in general change with time. In other words, the boundary dO,t moves. It is often 
known as the moving-boundary or free-boundary. 

Our work is motivated by the work of Caffarelli and Vazquez [2j in which authors studied 
this problem in the case p = 2. Their result stated essentially that if d^o € C^, /o € 
C'^{Qq) and A/q < 0, then there exists a solution to the problem. Moreover, if is 
compact, solutions vanish in finite time. Still in this case, long-time existence, uniqueness 
and regularity of the free-boundary have been studied by Daskalopoulos and Ki-Ahm Lee [5] 
or Petrosyan [II1II2] when the initial value is concave or star-shaped with bounded support. 
Other kinds of solution have also been studied (see also [9]). In the case p > 2, an elliptic 
version of the problem has been studied before by Danielli, Petrosyan and Shahgholian [3] 
or Henrot and Shahgholian [71|8]. As far as the parabolic problem when p > 2 is concerned, 
the only result we are aware of is by Akopyan and Shahgholian [1] where authors showed 
the uniqueness under the hypotheses that the time-section 0^ is convex and non-decreasing 
in time. The questions of existence or regularity of the free-boundary were not addressed 
in that paper. 

The main result of our work is stated below. 

Theorem 1.1. Assume that is a bounded and convex domain. The function /q is positive 
and concave in f^Q. Furthermore, on the boundary dO,Q, /o satisfies 

/o(x) = for all x 

l^/o(x)| = 1 for a.e. x. 

Then the problem ([P|) has a unique solution up to a finite time T where it vanishes identically 
in the sense that 

lim/(x,t) = VxGM". 
Moreover, the free-boundary dQt is smooth for all t G (0, T). 

It is well-known that solutions of the evolution p-Laplacian are only C^'" at points of 
vanishing gradient (see for example [6j). Hence, solutions to the problem ([P]) must be 
defined in some weak sense. We will state precisely the meaning of our solution in section 

El 

Our approach to the problem is totally different from [2]. To deal with the degeneracy, 
we will approximate the p-Laplacian with the following regularized operator 

(P(e)) A;f = dw{{\Df\' + ey-'Df). 

Here and throughout this work, we define q = p/2. We will establish some properties for 
solutions of these regularized problems and then let e go to to obtain a solution to the 
degenerate problem. 

In order to solve this regularized free-boundary problem, we employ a change of coordi- 
nates that transforms it into a quasilinear equation with Neumann boundary condition on 
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a fixed-domain problem. Applying results from standard theory of quasi-linear parabolic 
equations with oblique boundary condition, we show that this new problem admits a solu- 
tion for some positive time. Revert back to the original coordinates, we obtain a short-time 
existence result for the regularized problem. This argument is carried out in section [3l 

In section m we prove a simple estimate for the gradient \Df\ of solutions of the problem 



(P(e) ). In section O we prove a crucial result that the time-section fit remains convex and 
the function f{.,t) remains concave on Qt for all time t. Convexity of guarantees that 
the free-boundary dflt does not touch itself and also enables us to prove the non-degeneracy 
of 1-0/ 1 near the free-boundary. 

In section [6l we obtain an estimate for higher derivatives of / in a neighborhood the free- 
boundary d^lf, uniformly in time t and especially, in e, using the non-degeneracy of \Df\. 
This fact and the convexity guarantee that singular cannot develop on the free-boundary. 
The uniqueness for this regularized problem is obtained in section [71 In section [8l we then 
obtain a long-time existence result for solution of the regularized problem. Passing e to 0, 
we then obtain a solution to the degenerate problem in section [9j The uniqueness for the 
degenerate problem is then shown in section [TOl In the last section, we show that solution 
to our degenerate problem vanishes in finite time. 

Acknowledgement. I express my gratitude to my thesis advisor, P. Daskalopoulos, for 
suggesting this problem, and for her invaluable advices and support during the completion 
of this work. 



2. Definition of Solution 

In this section, we will define precisely what we mean by solution of the problem (P). We 
start by introducing some notations. For any < ti < t2 < T, define 

= nn{ti<t< t2}. 

First, we require that the free-boundary d^lt is in and the function / is in 

C{0,T;CHTk)). 

The equation 

ft = Apf in O 

is then defined in the sense that for any test function 6 in Cq°{Q) and for any < ti < t2 < 




\Df\P"'^ Df ■ DOd^dt. 



The Cauchy-Dirichlet conditions / = on d^t and /(.,0) = /o are understood in the 
pointwise sense 

/(x, t) — 5- as X ^ xq E dVLti 

/(x,t)^/o(x) ast^O. 
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Finally, the Neumann's boundary condition \Du\ = 1 is defined in the following classical 
sense 

fu{xo,t) = hm = 1. 

where xq is a point on the free-boundary dQt and i' is the spatial inward unit normal vector 
at Xq with regards to dClf 

3. Short-time Existence for Regularized Problem 

In this section, we will prove that the regularized free-boundary problem admits a solution 
for some positive time. We do it by a change of coordinates technique that transforms the 
problem into a fixed-domain problem. This technique has been used by other authors for 
different problems before (see for example [5], [4]). Note that concavity is not needed in 
this result. 

Lemma 3.1. Assume that C°° . The function /o is in C°°($7o) '^n-d positive in Qq. 

Furthermore, on the boundary O^Iq, /o satisfies 

/o = and \Dfo\ = l. 



Then there exists a smooth solution to the regularized problem (P(e)) for some T > 0. 



Proof. The argument in this proof works for any dimension, but due to the complexity of 
some computation involved, we will present the proof for the case n = 2 only. 

A word on notation used in this proof : we use bold-face letters x, y, ... to denote points 
in Euclidean spaces while normal letters x,y,z, ... for real numbers, scalars or components 
of points in Euclidean spaces. 

Denote by S the smooth surface z = fo{x,y), {x,y) £ flo. Let T = {Ti,T2,T^) be a 
smooth vector field on Jig such that T{x, y) is not a tangential vector to the surface S at 
the point /o(2;, y). Since |-D/o| = 1 on the boundary BVIq, we can also choose T to be parallel 
to the plane z = in a small neighborhood of OVLq. 

It is known that for some positive, small enough r/, we can define a change of spatial 
coordinates 

^> : X [-7/, rf\ 

by the formula 

The map <^ defines x, y and z as smooth functions of n, v and w with smooth inverses. 

The graph of {x,y, f{x,y,t)), {x,y) £ is then transformed to {u,v, g{u,v,t)), {u,v) G 
JIq via this coordinates change for some uniquely-defined g if the surface z = f{x,y,t) is 
sufficiently close to S ((x, y, f{x, y, t)) € ^{^o x [—r], rj]) for all (x, y) G Qt)- When / evolves 
as a function of {x,y), g evolves as a function of {u,v). Importantly, the domain of g is 
fixed as due to our requirement that T is parallel to the plane z = on OQq. 



x^ 






y 




V 


X ) 




[w 
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We will compute the evolution equation and the boundary condition of g. Denote by rc„, 
Xv, Xw, yu, Uv, Uwi Zu, Zy and the partial derivatives of the functions x{u,v,w), y{u,v,w) 
and z{u,v,w). Similarly we denote partial second derivatives of x,y and z by XumX^v, ■■■■ 

We begin with first derivatives. Since x, y and z are functions of u, v and w, while 
w = g{u, V, t) is a function of u, v and t, we have 

'dx dy_^ 
du du 

dx dy 
^dv dv ' 



Xu 


I ™ dw 
-1- Xu, 


Vu ~\~ Uw gy^ 


Xv 


I ™ dw 
-T Xw g^ 


Uv + Vw'g^ 


Xu 


~\~ X-wQu 


Vu + ywgu\ 


Xv 


~l~ XyjQy 


Vv + VwQv / 



We can compute the partial derivatives of u{x,y,t) and v{x,y,t) by 



(3.1) 



du 
dx 


du\ 
dy 1 




/dx 
1 du 


dx 
dv 


dv 
"Sx 


dv 1 


\-\ 


i dy 
\du 


dy 
dv 



-1 



dy 
1 / dv 



D 
1 

D 



_dy 
du 

Vv + VwQv 

~yu Uwdu 




{XuVv X^yu) ~\~ {Xuiyv UwXv^du ~l~ {yvjXu X^jyu^Qv- 



' du dv 

dx dx 

du dv 

.dy dy 





Next we compute the second-order derivatives. First, we have partial second order deriva- 
tives of x with regards to u and v. 



d'^x dw / dw \ ^ d"^ 

XuU ~1~ ^Xuyjj— -\- Xintn I 1 ~1~ X 



W 



-^UU I '^•^UW r^ I ^WW I o / ' UJ /-^o 

a^u ou \ou J a^z 

— Xuu ~\~ '^Xu;u9u ~l~ XyjQu 



d'^x 
d'^x 



Xyv ~1~ '^Xi^vQv ~t~ XyjQ', 



vv 



XuV XyjuQv XyjuQu "t" X^jfjuV 



dudv 

and similar formulae for y and z. 

Differentiate (j3.3p we have 

^ d'^z dz d'^u ^ dz d'^v ^ d'^z / du\^ ^ ^ d'^z du dv ^ d'^z / dv\' 
d'^x du d'^x dv d'^x d'^u \dx J dudv dx dx d'^v \dx J 
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We need to compute second order derivatives of u and v with regards to x and y. The 
formula (|3.4p is true if we substitute any function of u and v in place of z. Because second 
order derivatives of x and y with regards to x are zero 



^ dx d'^u dx d'^v d'^x 

du d'^x dv d'^x d'^u 

^ dy d'^u dy d'^v d'^y 

du d'^x dv d'^x d'^u 



du 
dx 



du\' 
dx J 



2 d'^x du dv d'^x f dv 
dudv dx dx d-^v \dx 



^ d^y du dv d^y 
dudv dx dx d'^v 



dv\' 
dx J 



In other words 



dx dx \ 




fd'^u 






d'^x 


dy dy 


II 


Xd'^v 


du dv ' 




^ d^x 



+ 



/ d'^x (du\'^ , 

d'^u \ dx) ' dudv dx dx 



du dv I d^ 

d'^v \dx) 



(dv-s'^X 



, d'^y /au\2 ^ d'^y du dv , d'^y ( dv\'^ 
\Wu \dx) '^dudv dx dx "i" Wv \dx) / 



or 



d'^u\ 




/dx 


dx 


d'^x 




du 


dv 


d^v I 




, dy 


dy 






\du 


dv 






/ du 


du 






1 dx 


dy 






[ dv 


dv 






\dx 


dy 



dx dx\ ^ / d'^x ( du\'^ , 2-^^^ — — + '^^^ ( dv\'^ 
du dv \ d-^u \ dx ) dudv dx dx d-'v \ dx / 

d'^y (du\'^ _|_ 2 d^y du dv , d'^y ( dv\'^ 
\ d'^u \ dx) dudv dx dx d^'v \dxi 



I d'^y (duy _|_ 2 du dv , d^y ( dvy J 

\ d-'u \dx) dudv dx dx d-^v \dx ) / 



d'^y du dv_ j_ d'^y ( dvy 



We then have 



dz d"^! 
du 3^3 



^ dz d'^v 
dv d'^x 



(dz 
\du 



( dz_ dz 
V du dv 



{fx fy) 





( d^x (du_~^ _j_ 2 Si^^ du dv , d^x f dv^ 

d-'u \ dx ) dudv dx dx d-^v \ dx ) 

d'^y ldu\'^ _|_ 2 d'^y du dv , d^y r dvy 

\ ff-'u \dx) dudv dx dx d-^v \dx ) 



( d'^x (du\'^ I 2J 

d'^u \dx ) dudv dx dx 



du dv I 



d'^x (dv\'^\ 
d'^v \dx) 



\ d'^y (du\'^ _|_ 2 d^y du dv , 
\ d'^u \ dx ) dudv dx dx 



d^y 
d-'v 



dvy I 

dx) / 



Substitute into (IHI 



fx 



d'^z ^ d'^x 
d'^u d'^u 



d'^y 



+ 2 



\dudv 



fyQ2^ 1 ( o„ 



du 
dx 



+ 



d'^z ^ d'^x 



d^y\ f dv\' 
dx J 



fyQ2y 



^ dudv ^ dudv J dx dx' 
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Let 



— fx^w fyV 



w 



u 



— Effuu ~l" '^{Z'oju fxXwu fyywu)9u "I" {^uu fx^uufyUv 

d'^z d'^x d'^y 
-D = 7^9 Jx-^ Jy'^ 



g2y Q2y y Q2 



V 



— ^9vv "I" ^{^wv fx^wv fyywv)9v "I" {,^vv fx^vvfyUvv) 

^ ~ o ci ~ Ja; o_ ci_ ^ 



then 



dudv dudv dudv 

EQuv ~I" (^«it) fxX'ojv ~ fyywv)9u "I" (•2-u;u fx^wu fyywu)9v 
(-^uii fx^uv fyyuv)j 



\dx J \dx J dx dx 

^({duV (dvV dudv \ F 

= ^ \ \ ] 9uu+ \ TT ] 9vv+ ^Tr^9uv + 7^ 



ax / \dx J dxdx ID'' 

where F is a smooth function of u,v, 9, Qu, 9w We have similar formulae for fxy and fyy 



du\'^ f dv\'^ ^dudv \ F 



_ / dudu dv dv du dv dv du\ \ F 

fxy — E I -pr--pr-9uu + 'K~1^9vv + 1 o ) 9uv J + 

\dx dy dx dy \ dx dy dx dy J J 

where F denotes different smooth functions of {u,v, g, Qu, 9v)- 
To compute ft, we differentiate z = f{x,y,t) 

dw 

Zw-Q^ = ft + {fxXw + fyyw)wt 
ft = (Zw - fxXw - fyyw)9t = Egt. 

Substituting into the equation for / 

ft = {\Df\^ + ey-^Af + {p- 2)(|i?/p + ey-\fxxf! + fyyfl + 2/.,/.^) 
and simplifying E from both sides we then obtain an evolution equation for g in the form 

gt = A'^{u,v,g,Dg)gij + B{u,v, g, Dg). 
On the other hand, the boundary condition \Df\ = 1 becomes 

C{u,v,g,Dg) = 

for some function C. 
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We claim the following is true when g = (i.e at t = 0) : 

• A^^ ,B and C are smooth functions of u,v,g and Dg. 

• (A^^) is positive definite. 

• C is oblique. 

Because the surface S and the vector field T are both smooth, it is clear that A^^ , B and 
C are smooth functions of u,v,g and Dg whenever 

„ dx dy dx dy , 
D = / 

du dv dv du 



E 



fxXw fyVw 7^ 0' 



The condition that E ^ Q follows from our choice that T is transverse to S. The condition 
D 7^ is a consequence of the fact that the function <I> is invertible in a neighborhood of S. 

Next, to show that {A^^) is positive definite, we write 

A'^ = {\Df\^ + tf-^Al + {p- 2){\Dff + e)''-^A'^ 

where A^ is the coefficient of gij (i, j G {u, v}) obtained from the transformation of A/ and 
A'^ from f^f^^ + f^fyy + 2fj;fyf^y. Wc cau compute exphcitly 

2 



A, 



+ 



dy) 



du dv I du dv ^ 
dx dx dx dy 



du dv I du dv_ 
dx dx dx dy 



'dv\'^ 
k dx 



Ao 



f Oil _i_ f du\ ' 
J^dx ~^ Jydy) 



( f J- f ^n. \ ( f ^ -I- f ^ 

WJ^dx ~^ JVdy ) y^dx JVdy 



f 9u I f du\ ( f dv_ I f dv\\ 
dx ~^ Jy dy I xJ^dx ~^ Jy dii I 

f §IE J- f d"" 1 

dx ~^ jy du I 



' dx J y dy j 
2 



'yay) 



It is obvious that Ai and A2 are non-negative definite. Furthermore, if D 7^ 0, Ai is actually 
positive definite (det(74i) = D^). It then follows readily that {A^^) is positive definite. 

For the proof that C is oblique, we refer to the Appendix of [5]. 

From the continuity, there must exist a positive number 5 such that those three claims 
are true for all g that satisfies |5'|(7i{no) < ^- then a consequence of standard theory 

of quasilinear parabolic equation with oblique boundary condition (see for examples jlOj . 
Chapter 14) that there exists a solution (7 up to a positive time T to the problem. 

gt = A^gij+B inl7ox(0,r) 
C{u,v,g,Dg) ={) on9J7ox(0,r) 
[ff(.,0)=0. 

This solution is actually smooth up to the boundary for all t G [0, T) since VLq is smooth 
and C{u,v, g, Dg) is a smooth function of {u,v, g, Dg). Choose a number T' in (0, T] such 
that \g\ < 1] on X (0, T'). Reverting back to the original coordinates system we then 
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obtain a solution to the regularized problem (P(e) ) up to T' . It is clear that the domain Q. 



is smooth and the solution / is smooth up to the free-boundary for all time < t < T' . □ 

4. Gradient Estimate 
Lemma 4.1. Assume the same hypotheses as in the Lemma \3.1[ Furthermore, /o satisfies 



\Dfo\ < 1 on rio- If f is a solution of the problem (P(e)), then 

\Df{x,t)\<l 

for all (x, t) e fl. 

Proof. We will show an equivalent fact that 

/A(x,t) < 1 

for any unit vector A. 
Let 

a'^iDf) = {\Df\^ + + 2(g - + e^-'^fj, 

where 6ij is the Kronecker delta function. Recall that we define q = p/2 throughout this 
work. Then the evolution equation of / can be written in non-divergent form as 

ft = i\Df\^ + ey-'Af + 2(g - l)i\Df\^ + ef-^f^.f^j 
= a''f^r 

We compute the evolution equation for f\ 

fxt = a'^fx^j + {fiJDa'^)Dfx. 

Since this equation satisfies the Strong Maximum Principle, f\ must attain its maximum 
value on the parabolic boundary of Q. Because /a < 1 on the parabolic boundary of 0,, it 
then follows that 

/a<1 

in 0, for all unit vector A. □ 

Lemma 4.2. Assume the same as in the last lemma, then at any point xq on the free- 
boundary dO-t 

fuuixo,t) < 

where v is the inward normal vector at xq with regards to dQt ■ 

Proof. Apply Hopf 's Lemma to the evolution equation for fi, from the last lemma, observing 
that fu attains the maximum value of 1 at (xo,t). □ 
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5. Convexity 

In this section we will show that the time-section Qt remains convex and the function 
f{.,t) remains concave on Qf Normally, for this kind of question, the main difficulty lies 
in showing that remains convex. The arguments for the case |) = 2 as in [11 j or [5j do 
not translate directly to the case p > 2. On the other hand, our argument here can be 
simplified to give a new and simple proof for the case p = 2. The argument relies heavily 
on the Neumann boundary condition \Df\ = 1. 

Lemma 5.1. Assume the same hypotheses as in the Lemma \3.1[ Furthermore, assume 
that 0,Q is strictly convex and /o is strictly concave on ^Iq. If f is a solution to the problem 



(P(e)) up to some positive time T, then is strictly convex and f{-,t) is strictly concave 



for all t G [0,r). 

Proof. We will show that 

/AA(x,t) <0 

for any point (x, t) G ^, any unit vector A and any t £ [0, T'] where T' is any number 
strictly less than T. Clearly this implies that fit is strictly convex and / is strictly concave 
for ah t G [0,T). 

First, we compute the evolution equation of /a a, 

ft = a^'fij 

fxt = a'^fxij + ifijDa'^) ■ Dfx 

fxxt = a'^fxx^j + 2ifx^jDa'^) ■ Dfx + ifijDa'^) ■ Dfxx + {fij{Da'^)x) ■ Dfx 

= a'^fxMj + {f^JDa'^) ■ Dfxx + {2fxijDa'^ + fijD{a'^)x) ■ Dfx 

Since / is smooth for all t G (0, T), there exists a finite number C{T') such that for any 
unit vector A and any point (x, t) G H {0 < t < T'} 

\2fxijDa'^ + fijD{a'')x\<C. 
Choose a smooth function vq on M" such that 

fo < < -(/o)aa in ^^o 
[ vo = on dQ.Q. 

Such vq exists because /o is strictly concave on VLq. Let v be the solution of the Cauchy- 
Dirichlet problem 

vt = a^^Vij + (fijDa^^) ■ Dv — Cv in Q 
v = ondnn{0<t<T} 
v{., 0) = vq in r^o- 

Applying Strong Maximum Principle and Hopf 's Lemma to v we easily deduce that 

{V > m Q 

\Dv\>0 ondnn{0<t<T}. 
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We are going to show that 

(5.1) v + fxx<0 

for ah t G [0, T'] and all unit vector A. Assuming that it is not the case, i.e there exists 
some point (x',t') and some unit vector A' such that 

(t; + /vy)(x',O = 

and 

v + fxx<0 

for all t < t' and all unit vector A. In other words, t' is the first time ()5.ip fails. We consider 
two cases, (x',t') is an interior point or a boundary point. But first, note that we have the 
evolution equation for V = v + fyx' 

(5.2) Vt = a'JV-, + (^Da'J) • DV + {2fxHjDa'^ + hD{^')x') ■ Dfx' - Cv. 

If (x',f') is an interior point, then because it is a maximum point of V in ^It', we have 

a'^Vij < 
DV = 0. 

Substitute into <\5.2h we have 

Vt < i2fxujDa'^ + fijOia'^Yx) ■ Dfx - Cv. 
Because at the point (x',t') 

V + fx'X' >v + fxx 

or 

/a'A' > fxx 

for any other unit vector A, we have 

fxx' = 

for any A±A'. Hence, 

Vt < {i2fxijDa'^ + fijDia'')x) ■ A') fx'X' - Cv 
< —Cfx'x' ~ Cv (remember fyx' = —v < 0) 
= 

which contradicts the assumption that (x',t') is the first time V = 0. So (x',t') cannot be 
an interior point. 

If x' is on d^t'- Again, denote by the inward normal unit vector to d^t' x' . Then 
at this point we have from definition of (x',t') and A', 

{v + fyx')u<0 

fuX'X' < -Vu < 0. 

We will show that on the other hand 

(5.3) Ux'X' = 0. 
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We have from the Lemma 14.21 that 

fvv < 0. 

We also have as a consequence of the fact that \Df\ = 1 on the free-boundary and \Df\ < 1 
in the interior that fu\ = for any tangential unit vector A. Hence as a consequence of 
the fact /a'a' = 0, A' must be a tangential vector of dO,t'- Otherwise, there would be a 
tangential vector A that lies on the same plane with and A' such that 

/aa>0 

which contradicts our assumption on (x',t') and A'. 

Without loss of generality, we can assume that = ei and A' = 62- Because ei is the 
unit normal vector of dilt' at x' , in a small neighborhood of x', we can write dilt' as the 
graph of a smooth function 

xi = 'y{x2,x') 

where x' = (xs, ...Xn)- From here to the end of the proof, we will use 7' and 7" to denote 
the first and second derivatives of 7 with regards to X2 ■ Differentiate / = with regards to 
62 we have 

/17' + /2 = 

or 7' = since /i = 1, /2 = 0. Differentiate one more time and disregard all terms 
containing 7' we have 

fll" + /22 = 

and so 7" = since /22 = due to our assumption. Differentiate Df ■ Df = 1 twice with 
regards to 62 and disregard all terms containing 7' or 7" we obtain 

Df ■ Df22 + Df2 ■Df2 = Q 

or 

/l22 + |i^/2p = 0. 

As above, because fii<0 = ^22 ^i, we have 

f2i = 

for all i / 2. But /22 = as well, so D/2 = 0. Hence 

/122 = 

which is exactly what we want to show in (j5.3p . We then have a contradiction. In other 
words 

/aa<0 

for all {x,t) € $7 n {0 < t < T} and all unit vector A or equivalently, ^It is strictly convex 
and f{-,t) is strictly concave in ^It for all t € [0,T). □ 



a free-boundary problem for the evolution p-laplacian 
6. Regularity near the Free-Boundary 



13 



In this section, we show that the degeneracy \Df \ = is kept away from the free- 
boundary. Consequently, the free-boundary is smooth, uniformly in e. It enables us to show 
that the limiting function obtained by letting e go to satisfies the boundary condition of 
the original problem. The proof depends crucially on the concavity of /. 

We introduce some notations. We denote by i?r(x) the disk of radius r around x 

Br{x) = {y G M" I |y-x| < r} 

when X € and r E M. We write Br for Br{0). We also define 

Ar = {(x, t) I dist(x, dQt) <r}nn. 

For any point x = {xi,X2, ■■■Xn), we define 

V'i(x) = Xi 

^'(x) = {X2, ...,Xn). 

Lemma 6.1. Assume all hypotheses as in the Lemma l5.il Assume also that there exist 
positive numbers r, R and m and a point xq such that 

(6.1) 5r(xo) C J7t C 5r(xo) for all t G [0, T) 

(6.2) /(x,t) > m for all (x,t) G Br (xo) x [0,T). 

Then for any < Ti < T and k £ Z+, there exist positive numbers d{r, R,m,Ti) and 
C{d, k) such that 

Proof. To simplify the notation, we assume that the conditions (|6.ip holds for xq = 0. In 
other words 

Br cntCBR for all t G [0, T) 

/(x, t)>m for all (x, t) G Br x [0, T). 

Let {P,t) be a point on 9$7 for some t G [Ti,T). Fix this value of t from here until the 
end of this proof. Without loss of generality, we can assume that 

V'i(^') <o 

ip'{P) = 

First, we will show that /i(x, t) is bounded away from in a neighborhood of P in i}f 
Consider any point Q in ^If that satisfies the following conditions 

MQ)<o 

|V'(Q)| <r 
f{Q,t)<m/2. 



14 



TUNG TO 



Let R = (0, 'tp'{Q)). Because R G B^., we have f{R, t) > m. We also have /(Q, t) < m/2 and 
fi (x, t) decreases in xi as a consequence of concavity (here /i denotes the first derivative 
of / with regards to xi). Thus, 

m/2 < f{R,t) - f{Q,t) 

fi(,{xi,^l;'{Q)),t)dxi 



< \MQ)\fiiQ,t) 

<Rfi{Q,t) 
h{Q,t)>m/2R 

We just showed that if x satisfies 

^i(x) < 
IV''(x)| < r 
/(x,t) <m/2 

then 

fi{^,t)>m/2R. 

On the set containing all such x, the Implicit Function Theorem says that there exists a 
function g defined on the set 

B = {{y, x') €Rx I < y < m/2, \x'\ < r} x [0, T) 

such that 

^i(x)=5(/(x,f),V''(x)). 
We will compute explicitly the evolution equation and boundary condition of g 

91 

f - 

Ji — 

91 

f 9t 

It — 

91 

f 911 
111 = 3" 

9i 

f _ 9i9ii - 9i9ii 

Jli — 3 

91 

0? Qii — Q^ Qi Q^ i - Q^ Qi Oi + Qi Q-i Q^ i 
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The boundary condition \Df\ = 1 on dCl^ is equivalent to 

/ n \ 1/2 

on 

{{y,x') G R X R"-^ \ y = 0, \x'\ < r} X [0,T). 

Next we compute the evolution for g onB. In all ^ appearing in the following computations, 
unless explicitly marked otherwise, indices i and j run from 2 to n. Let 



Then 



fifjfij = — \ (911 -'^^9i{9i9ii - 9i9n) 
9i ^ 

+ X 9i9ji9i9ij - 9i9i9ij - 9i9j9u + 9i9j9n)) 
= — \ (giiM^ - 2giM X 9i9u + 5i X 9i9j9ij) ■ 



9i 

Substitute into the equation for /t, 



^511 + Pi X 3ii - 2^1 X 1^) 



9l 

+ 2(g - 1) + (^iiM^ - X + 5? X 

We want to show that there exist positive numbers A, A, independent of e and t such that 

A < V^iiij < A 
for all unit vector ^ G M". First, because 

1 > l^/l >fi> m/2R, 

we have 

1 < < £fi < 2i2/m. 

The upper bound A then is obvious. For the lower bound, because 
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it is enough to show that 

- 251 Yl 9iiiii +giY.^t>^ 

for some positive A. We have 

since g'l> M = \ + Y^gf and so 

^ + 5- ) il - 2gmCiCi + (gj - > for all 2 < i < n. 

Summing up we obtain 

- 291 Y, a,(i(i + 9? E a 

From the theory of quasi-linear parabolic equation with oblique boundary condition we 
can choose d < min(r, m/2) such that g{.,t) is in C°° on the set 

{{y,x') G M X M"-^ \ 0<y<d, \x'\ < d} x [Ti,T) 

and for any k, the norm \g{.,t)\Qk depends only on k, d,r, R,m and Ti, not on e, t or g^. 
Revert back to /, we conclude that / is smooth on the set 

{(x,t) G n I iv^'lx)! < d,TPi{x) < o,/(x,t) < d} n{ri < t < T} 

and again, for any k, the norm \ f{.,t)\Qk on this set depends only on k, d,r, R,m and Ti. 
Note that the above set includes the set 

B{P,d)rint. 

The conclusion is of course true for any point on dflt in place of P where t G [Ti,T). The 
lemma then follows. □ 



7. Comparison Principle 
In the first lemma here, we show that if /q is strictly greater than /q, then a solution to 



the problem (P(e) ) with initial value /q remains strictly greater than a solution with initial 
value /o. 

Lemma 7.1. Suppose that f and f are solutions up to some finite time T to the problem 



P(e) and P(e') respectively for some e > e' > 0. Suppose also that at the time t = 0, 

/o(x) < /o(x) in Qq. 
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Then for any t E (0, T), 

/(x,i) </'(x,t) inWt. 

Proof. Since 

/o(x) < /|^(x) in Ih 
we can choose a positive number m such that 

/o(x) + m < /o(x) in Qq. 
Choose a positive number 6 such that 6T < m. We will show that 

f'(x,t) - f{x,t) - m + 5t > inTTt 

for all t € [0, T). Assuming it is not the case, there must be a first time to such that at least 
one of the two above conditions is violated. Assume that the first condition is violated at 
to- In other words, dO,t and dQ[ touches at some point xo, then at that point 

/'(xq, to) - /(xq, to) -m + 5to = -m + 5to < 

which implies that the second condition must be violated before time to, contradicting our 
choice of to- Hence 

for all t G [0,to]. The second condition is violated implies that there is a point xq G i^to 
such that 

/'(xo, to) - /(xo, to) - m + 6to = 0. 

We consider the case xq S f?J7(g C first. Let u be the inward unit normal to dQto at 
Xo- From the definition of (xo,to) we must have 

(/'(xo, to) - /(xo, to) - m + 6to)^ > 

fu{^o,to) - ^(xo,to) > 
/^(xo,to) > 1 

which contradicts the result of Lemma 14.11 

If Xo G C f^to' thsn because it is an minimum point for f' — f on f^to) we have 

Z)/'(xo,to) = Z)/(xo,to) 
0> A/'(xo,to) > A/(xo,to) 
> /^^(xo,to) > /!.,.(xo,to). 
Plug into the equation for // and ft, recalling that e > e' we obtain 

/; = (|D/f + e'y-'^f + 2{q - l){\Df'\^ + e')"-'!/?/?/.. 
> {\Df\^ + ey-'Af + 2{q - l){\Dff + ey-^lDfll, 
= ft. 
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On the other hand, because is the first time f — f — m + 5t = 0, 

n-ft+s<o. 

Again, we arrive a contradiction. In other words, 

f -f -m + dt>0 

for all t € [0, T). The Lemma then follows readily. □ 
Remark 7.1. If we let 

m min {/o(x) - /o(x) | x G SIq } 

we actually prove that 

m{t) = min {/'(x, t) - /(x, t) | x G } 
is a non-decreasing function. 

We prove a slightly improved version of the last lemma. 

Lemma 7.2. Suppose f and f are solutions up to time T to the problem P(e) and P(e') 
respectively for some e > e' > 0. Suppose also that at the time t = 0, 

/o(x) < /o(x) in r^o- 

Then for any t G [0, T), 

Qt c n't 

f{^,t)<f{^,t) in^t. 

Proof. Without loss of generality, we can assume that /o attains its maximum value at the 
origin. For each positive A, define 

VL^ = {(x,t) I {\^,\h) G 0} 
/^ = ^/(Ax,A2t) 
1^0 = {x I Ax G ilo} 
/o^ = i/o(Ax,A2t). 



It is clear that / is a solution to the problem (P(e)) with respect to the initial data /g. 
Furthermore, since /o is concave, for each A > 1 we have 

/o^</o</o inf^^ 
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The last lemma says that for all t € [0, T/A^), 

^/(Ax,A2t) = /^(x,t) </'(x,t) in^^^ 

Let A ^ 1 we obtain 

Qt c n't 

/(x,t) </'(x,t) inQt 
for alH G [0,r). □ 



Corollary 7.3. Suppose (/, f^) and {f',0,') are two solutions to the problem (P(e)) with 
respect to the same initial value Jq up to time T. Then f = f on M" x [0,T). 



8. LoNG-TiME Existence for the Regularized Problem 

Lemma 8.1. Assume that /o satisfies all hypotheses of the Lemma \3.1\ Then there exists 
a unique solution to the problem P(e) up to some positive time T > where 



lim /(x,t) = 



for all X G 



Proof. Let T be the maximal existence time for solutions to the problem P(e) with the 
initial data /q. Due to the uniqueness result in section \T0\ there must be a solution / that 
exists up to time T. From the short-time existence result, T must be positive. We will show 
that 

lim/(x,t) = forallxGM". 

Assuming otherwise, then the same argument in the proof for the Lemma 111. II can be used 
to show that T must be finite. In other words, due to the concavity of /o, there exists a 
number c < such that 

div {{\Dfo\^ + ey-^Dfo) <c in no 

and consequently, 

^ ^ max/o 
~ |c| 

We will prove that we can then extend this solution to a time T' > T. From the concavity 
of f{.,t), / is a decreasing function in t. Define 

/t(x) = lim/(x,t). 

Since \Df \ < 1 in 17, fx is continuous. Because /r is not identically 0, there exist a ball 
Br{x') and a positive number m such that 

fT>fn in Br{x'). 
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From the Lemma |6.H for all t € [T/2,T), there exists a positive number d such that / is 
smooth up to the boundary and time T in the set 

{(x,t) I dist(x,9J]t) < d} nO[T/2,T)- 

Combine with the smoothness (depending on e) of / up to time T in the interior of ^[t/2,t) 
from the standard theory of parabolic equation, we obtain the smoothness up to the bound- 
ary and time T of / in ^[t/2,t)- Consequently, fx is smooth up to the boundary. From 
the Lemma 15.11 we know that ilx is convex and /t is concave in Qj-- However, we need 
a stronger result that is strictly convex and fx is strictly concave in in order to 
apply the Lemma l3.1i In deed, we can improve the result in the lemma [5TT] bv duplicating 
the proof and substituting T' by T directly. In that proof, because we did not have the 
smoothness of / up to time T, we need to introduce T' < T to guarantee the existence of 
a finite number C{T') such that 

\2fxijDa'^ + fijD{a'^)x\ < C 

for all t S [t,T']. But now we have the smoothness of / up to time T, we can derive the 
fact that there exists a number C(T) such that the above inequality holds for all t € [0,T). 
The proof then guarantees that / is strictly concave at the time T. 

The function fx now satisfies all hypotheses of the Lemma l3. II By that Lemma, we can 
then extend the solution / to some time T' > T. It contradicts the maximality of T. So we 
must have 

lim /(x, t) = 

for all X e M". □ 



9. Existence of Solution to the ^-Laplacian problem 

In this section, we will pass e to and obtain a solution to our degenerate problem. 

Lemma 9.1. Assume that VLq is a bounded and convex domain. The function fo is positive 
and concave in Furthermore, on the boundary dO,Q, /o satisfies 

/o(x) = for all x 
\Dfo{^)\ = 1 for a.e. x. 

Then there exists a solution to the problem (jP]) up to some time T where 

lim/(x,t) = VxGR'". 

The free-boundary d^lt is smooth for all t € (0,T). 



A FREE-BOUNDARY PROBLEM FOR THE EVOLUTION p-LAPLACIAN 21 

Proof. Choose a sequence of functions /g with positive sets ilg for all e G (0, 1) such that 

G C7°° and G C°°(ll), 
rig is strictly convex, 
/o is strictly concave, 

C 17^^ and /o^i < /o^2 if ei > £2, 
J^o = UJ7^ and /o(x) = lim /o"(x) for all x G R", 

e— >0 

{Df'l = 1 on dft'. 

In other words, /q is an increasing sequence of smooth and strictly concave function that 
converge to /o as e ^ 0. From the Lemma |8.H for each e, there exists a unique solution 
to the problem (P(e) ) up to some time where it vanishes identically. We will prove that 
lim^^o is a solution to the original problem (jP]) . 

From the lemma [7^2] and our choice of /g, it is clear that if ei > £2, then 

f'^ < in 

Define 

T = lim 

n = un' 

/(x, t) = lim /'(x, t) for (x, t) G M" x [0, T). 
Due to the uniform smoothness of in a neighborhood of d^l^, we have 

/ = and \Df\ = 1 on 91) x {0 < t < T}. 

If (x, t) G ri, there exists an eo such that (x, t) G f^"^ for all e < eo- Since /"^(x, t) increases 
as e decreases to 0, 

/(x,t) =limr(x,t) >0. 

From the bound \Df \ < 1 and the Corollary 2.15 in Chapter II of [TC^, we have interior 

C^'^^'^ estim 
for every x, 



C^'^^'^ estimate for as functions of t, uniformly in t and e. Together with the fact that 



we have 



Because for every e and t 



lim r(x,t) = 



lim /(x,t) = 0. 



r(x,t) <r(x,0) </o(x,0). 
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we have 

/(x,t) = limr(x,t)</o(x). 
On the other hand, since /'^(x, t) increases as e decreases to 0, 

for any e. Consequently, 



lim/(x,t)>limr(x,i) = /o^(x) 



Hence 



hni/(x,t)>lim/o^(x) = /o(x). 



lim/(x,t) = /o(x). 



From \Df^\ < 1, we can choose a sequence of e converging to such that 

Dr ^ Df 

in all compact subsets of Jl. Given any function 9 € C^{fl) and < ti < t2 < T we have 
from the equation 

f,^ = div{{\Dr\^ + er-'Dr) 

that 

J f'6tdxdt- J f'6d:ic '= J {\D + ef-^ D f ■ DO d^dt. 

Passing to the limit we then obtain 

j fOtdxdt- J fOdx '= j \Df\^^''-^^ Df ■ DOd^dt. 

□ 

10. Uniqueness 

Lemma 10.1. Solution obtained in the Lemma \9.1\ is the unique solution to the problem 
(P)- 

Proof. Assume that there exists another solution g to the problem (P). Let J7* be the 
positive set of g and T* its existence time. Also without loss of generality, assuming that 
/o attains it maximum value at 0. For each positive A, it is clear that 

5"(x,t) = A-VA'x,AP+2t) 



is a solution to the problem P(e) with positive set 



= {(x,t) I (A^x.A^+^t) e 17*} 

and initial data 

go"(x,t) = A-Vo(A2x,A^'+2t). 
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Clearly, for A < 1, 

/o < 9o in 1^0- 
We will show that for ah t < min(r, X-(p+'^)t*), 

/(x,t) <5^(x,t) inQf 

Assuming it is not the case, then there must be a first time to where at least one of those 
two inequalities is violated. If the first one is violated at the time to, it means dflt and dQ^ 
touch at some point xq. At that point (xo,to)) 

|Z)/| = 1>A = |I)5^|. 

There must be then a point x G ilj such that 

/(x,to) > /(x,to)- 

which implies that the second inequality must be violated at some time before to. So, up 
to time to) 

Wt C Qt- 

Consequently, on the parabolic boundary of f^[o,to]i f < 9^- Thus, from the lemma 3.1 in 
Chapter VI of [6], we have f < g in ^[o,to] which contradicts our choice of to- Hence for all 
t < min(T,A-(P+2)T*), 

/(x, t) < g\x, t) = A-ig(A2x, Xf+H) in 
Let A — > 1 we obtain for all t < min(T, T*), 

/(x,t) < 5(x,t) in Of. 
Arguing similarly for A > 1 we on the other hand obtain 

ff(x,t) </(x,t) inQl 

Thus for all t < min(T, T*) 

/(x,t) =g(x,t). 



□ 
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11. Vanishing in finite time 
Lemma 11.1. The existence time T of the solution obtained in the Lemma \9.1\ is finite. 

Proof. Clearly from the Comparison Principle in section [7] and scaling that it is enough to 
prove this lemma for one particular initial function /q. We show that if a smooth function 
/o satisfies all hypotheses of the Lemma l9.ll and 

Ap/o < c 

for some c < 0, then for any < ti < t2 < ^ and any x E r^tj, the corresponding solution 
/ satisfies the inequality 

(11.1) /(X,t2)-/(x,tl) <c(t2-tl). 

It then readily follows that 

^ ^ max/o 
~ |c| 

Choose the sequence {/g} so that 

div((|D/o^|2+eriz)/o^)<c. 

We will show that satisfies 

rt<c 

for all e and (jll.ip then follows immediately. 

Differentiating the equation satisfied by with respect to t, it is easy to see that f^ 
satisfies the Maximum Principle. Since f^ < c at the time t = from our choice of /q, all 
we need to show is that f^ cannot attain its maximum value on the free-boundary. Prom 
the Hopf's Lemma, if f^ attains its maximum value at a point xq on the free-boundary 
dVtto^ then we must have 

(/^)t(xo,to) = (/n.(xo,to) <0 

where v is the inward unit normal at xq with respect to d^to ■ On the other hand, since fit 
shrinks in time, for any t < to G J^t and so, 

/^(xo,i) < 1 

while 

fui^O,to) = 1 

which lead to 

(/^)t(xo,fo) >0. 

□ 

Remark 11.1. The finiteness for the existence time holds for any initial data with bounded 
support, not just concave ones. 
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